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TECHNICAL NOTE 2817 


A THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE 
EFFECTS OF YAW ON PRESSURES, FORCES, AND MOMENTS 
DURING SEAPLANE LANDINGS AND PLANING 
By Robert F. Smiley 

SUMMARY 


A theoretical Investigation was made of the hydrodynamic forces and 
moments experienced diiring yawed water landings and planing of seaplanes 
of arbitraiy constant cross section. Equations are developed for the 
side force, the rolling moment, the yawing moment, the pressure distri- 
bution, and the peak pressure. For the special case of the non-chine- 
immersed straight- sided wedge, these equations are such that the time 
histories of the side force and rolling and yawing moments can be 
expressed as families of generalized curves. 

Experimental measiirements of the side force, the rolling and 
yawing moments, the pressure distribution, and the peak pressure, 
obtained d\iring landing and planing tests made in the Langley impact 
basin with a float having an angle of dead rise of 22.5°^ are presented 
and compared with the theoretical predictions. The landing tests cov- 
ered yaw angles between 0° and 12° for trims of 3.2°, 6.3°, and 9.3° 
and the planing tests covered yaw angles of 60 and 9° for trims of 6.3° 
and 9 • 3 °' In general, the experimental data appear to be in reasonable 
agreement with the corresponding theoretical predictions. 


INTRODUCTION 


In recent years much experimental and theoretical research has been 
directed toward the obtaining of information concerning the motions and 
hydrodynamic forces experienced during the landing or planing of sea- 
planes. For the case of the symmetrical landing or planing many theo- 
retical and experimental data are now available (refs. 1 to 22 ). For 
the case of unsymmetrical landing or planing, however, very little exper- 
imental. or analytical work has been done ( refs . 23 to 25) . 

The purpose of this paper is to approach this problem of unsymmetrical 
loads by a theoretical and experimental study of the effects of yaw on 
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seaplane loads and motions. The analysis deals, in general, with the 
treatment of yaw forces on hulls of arbitrary constant cross section 
and deals, in particular, with the case of a wide straight- sided wedge. 

A study is made of the pressure distribution, side force, and 
rolling moment occurring during the two-dimensional yawed Impact of a 
symmetrical body on a smooth water surface. The problem of the oblique 
yawed impact and planing of a three-dimensional body of arbitrary con- 
stant cross section is then solved by the use of the results of the two- 
dimensional analysis. After a consideration of several constants needed 
for theoretical computations, experimental impact and planing data 
obtained from tests in the Langley impact basin, which are described 
in the appendix of this paper, are compared with the corresponding theo- 
retical predictions. 


SYMBOLS 


Any consistent set of iinlts may be used unless otherwise indicated. 
Bars over symbols indicate the average value of the quantities involved. 


A hydrodynamic aspect ratio, 

(Wetted length at keel)^ 

Wetted area projected normal to keel 

perpendicular distance between keel and sixis of rolling 
moments 




distance between step and axis of yawing moments measured 
parallel to keel 


B, B]l, B2 side-force factors 


4 

I 


b 

c 


E, El, Eg 
Fi, F2 


beam of hull 

wetted semiwidth of hull in any transverse section 
rolling-moment factors 

total hydrodynamic force on two sides of a wedge 

hydrodynamic force normal to keel (parallel to ^-axis) 

hydrodynamic force perpendicular to seaplane plane of 
symmetry (parallel to T^-axis) 
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f ( P) dead-rise function 

G yawing-moment factor 

g acceleration due to gravity,- 32.2 ft/sec^ 

H perpendicular distance “between keel and plane of chines 


h 

J 

K 


L 

Mg 




ni, U2 


wetted height of hull in any transverse section 
empirical fimction of angle of dead rise 
constant for pressure calculations, 

K % i/l - 3 tan^P cos 3 _ tan p sln^3\ 

" 1 . 7^2 " 3 . 3 „ J 

longitudinal-wave-rise ratio 

hydro(3ynamic rolling moment about an axis parallel to l-axis 
and in -plane; the moment is considered positive if it 
■bends to decrease the angle of dead rise on side of hull 
having the largest Impact force (advancing side of hull) 

hydrodynamic yawing moment about an axis parallel to ^-axis 
and in ^5 -plane; the moment is considered positive if 
it tends to reduce angle of yaw 

two sides of a thick rectangular flat pla-te 


ni_ side-force load factor, ^ 

I w 

•• 

ni„ vertical load factor, ~ 

“ g 


p Instantaneous pressure 

R transverse-wave -rise ratio 


Such wetted area of one side of seaplane projected normal to 

plane of symmetry of seaplane at chine immersion 

« 

increase of wetted area of one side of seaplane projected 
normal -to plane of symmetry of seaplane subsequent -to 
chine immersion 
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S2 wetted areas of two sides of seaplane projected normal to 
the X-axis 

s distance between a given flow pl^e and keel-water- surface 

intersection (parallel to |-axis) 

t time after water contact 


h 


V 

W 

w 

X,Y,Z 


X',Y',Z' 


X 

z 

z 


• ■ 

z 

a 


3 


velocity of hull relative to undisturbed water 
velocity of hull normal to keel (parallel to ^-axis) 

acceleration of hull nomal to keel (parallel to ^-axis) 

velocity of hull perpendicular to seaplane plane of symmetry 
(parallel to T]-axis) . . 

local fluid velocity relative to body at body surface 

weight of hull 

thickness of a rectangular flat plate 

coordinate axes with respect to plane of seaplane landing; 
X-axis is in plane of seaplane motion and is parallel 
to water surface, Z-axis is in plane of seaplane motion 
and is perpendicular to water surface, and Y-axis is per- 
pendicular to X- and Z-axes 

X-, Y- , Z-coordinate axes rotated about the Z-axis through 
yaw angle; X'-axis is in plane of symmetry of seaplane and 
is parallel to water surface, Z*-axis is in plane of sym- 
metry of seaplane and is perpendicular to water surface, 
and Y'-axis is perpendicular to X’- and Z'-axes 

horizontal velocity of hull (parallel to X-axis) 

horizontal acceleration of hull (parallel to X-axis) 

vertical displacement of step (parallel to Z-axis) 

vertical velocity of hull (parallel to Z-axis) 

vertical acceleration of hull (parallel to Z-axis) 

3 B 

constant for free- lateral- mot ion solution', 

[f(p)] cp(A) 

angle of dead rise, radietns 


Pi 


1 PH 

average angle of dead rise, tan — , 

b 

cramma function 




S 


A 


radians 
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7 

€ 

• « 

9 

X 




I 

p 

T 

cp(A) 


fli^t-path angle 

transverse slope of hull surface at any point 
penetration of hull nomial to keel (parallel to ^-axis) 
approximation for V^, x sin t + i cos t 

approximation for x sin t + z cos t 

absolute value of transverse distance from center of 
model (parallel to T)-axis) 

effective angle of dead rise 

parameter used in equation ( 29 ) 

coordinate axes with respect to seaplanej the 5-axis is 
in plane of symmetry of seaplane and is parallel to keel, 
the Ti-axis is perpendicular to plane of symmetry of sea- 
plane, and the 5-axis is perpendicular to I- and T|-axes 

distance forward of step (parallel to 5-axis) 

mass density of water, 1-938 slugs/cu ft for tests 

trim 

velocity potential on surface of body 
aspect-ratio correction 
yaw angle 


Subscripts: 

a 

ch 

cp 

f 

max 

o 


restricted side motion 
at chine Immersion 
center of pressure 
free side motion 
maximum 

at water contact 





6 


NACA TN 2817 


P 

u 

1 

2 


peak 
unyawed 
normal flow 
transverse flow 


Dimensionless variables: 


Cd 


Ct 


sin T - . . 

approach parameter, cos 

vertical- load- factor coefficient, 
nlw8 Tgw sin T cos^t\ 

* 0 ^ 1 [f ( 3)3 ^cp(A) rtpg 


draft coefficient, z. 


[f(3^ ^(A)jtpg' \^'^^ 

sin T cos^ J 


time coefficient, tzo c ^ r 

1 6W sin T cos^ 


Ct^ or 


side-force coefficient (C^j for restricted case, for 

free case^, 


Il 


Birpxoto sin. i|r 




2 v7 tan t 


6w 


[f( 3)3 %(a) jtpg 


2/3 


Cm I 


rolling-moment coefficient for restricted cape, 
Mt 


Tr2BEff^Q sin f 

i;|f(3)32(p(A)g sin^3 


rolling-moment parameter, a 


12 sin^3 r rf(P)]^tp(A)n:pg 


jt^E 


6W tan T 


1/3 
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Cm^ yawing-moment coefficient for restricted case, 

sin t 

jf(3)[|^cp(A)g tan T 




yawing-moment parameter. 


„ 3 f [f(p)]%(A)3rgp tan^T 



L 6W . 


1/3 


ANALYSIS OF TWO-DIMENSIONAL PROBLEM OF YAWED IMPACT 


General analysis .- Before the specific problem of the yawed impact 
of a seaplane is discussed it is desirable to consider first the corre- 
sponding two-dimensional problem, namely, the oblique impact of a sym- 
metrical body on a smooth water surface (fig. l) . For convenience, the 
body is assumed to be stationary and the fluid is assumed to be moving 
toward the body at a velocity V far from the body, the vertical and 
side components of , this velocity being designated as and Vt|, 

respectively. Buoyant and viscous forces are neglected. The velocity 
potential for the flow and the corresponding fluid velocity at the body 
are designated as 0 and v, respectively. As a first crude approxi- 
mation this flow field can be considered to be the sum of the following 
two flows (see fig. 2): (l) the flow for the symmetrical vertical impact 

of the same body on a fluid moving at a velocity at infinity, which 

will be designated by the subscript 1, and (2) the symmetrical horizontal 
flow about the corresponding vertically symmetrical completely submerged 
body at a velocity which will be designated by the subscript 2. It 

should be noted that tjie application of this superimposition principle 
is not strictly correct because of the presence of the free surface in 
the basic problem. For lack of a better or simpler way to approach this 
problem, however, the principle is assmed to be approximately valid. 

According to Bernoulli's equation the pressure (above atmospheric) 
on the body for the three problems represented in figures 1 and 2 is 

p = ^p(v2 _ v2) - p H ( 1 ) 

PI = ip(vj2 . vp2) - p 


( 2 ) 
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P2 





(3) 


where equation (l) refers to the complete yawed flow of figure 1 and 
equations (2) and ( 3 ) refer to the vertical and horizontal components 
of the flow in figure 2, respectively. As the flow field of figure 1 
is assumed to he the resultant of the flow fields in figure 2, then on 
the bottom of the body 


^ 

V = V]^ + vg (5) 


The plus sign in equation (5) refers to the left side of the body and 
the minus sign refers to the right side (see arrows designating fluid 
flow directions in fig. 2). Substituting equations (U) and (5) and the 

relation into eqiiation (l) gives 


P = + V.,2 . . p ^ . P ^ 

= |p(v5® V,,2 - t - YgS) ■ ^ 


This analysis is now restricted to the determination of the first- 
order effects of the yaw; that is, the term in equation (6) which, is 
proportional to V|^ (that is, pvjV 2 ) is considered but terms proper- 


that is, ^Vt|^, and p 

Then equation (6) reduces to the relation 


are neglected. 


P 




PV1V2 


(7) 


The difference in pressure between that for the yawed case and that for 
the unyawed case can be obtained by subtracting equation (2) from equa- 
tion ( 7 ) and is 


P - Pi = i PV1V2 


( 8 ) 
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Equation ( 8 ) gives a first approximation for the pressure due to yaw on 
an arhitrarily shaped symmetric body. Because of the preceding simpli- 
fying assumption that a super imposition procedure is valid and because 
of the neglect of second-order terms proportional to the yaw velocity 
Vt), this equation should not necessarily be considered to be applicable 
to practical yaw conditions. Rather it is or shoxild be considered as 
an equation which is probably approximately valid for infinitesimal yaw 
angles and which may be corrected later empirically to apply to the 
practical case. 

According to equation ( 8 ) (see statement after eq. (5)) the unsym- 
metrical pressure or the pressure due to yaw is positive on one side of 
the body and negative on the other so that the net upward force due to 
yaw is zero. The vertical force on a yawed body at a vertical velocity 
thus is to a first-order approximation the same as that on an unyawed 

body at the same vertical velocity or 

F5 = ( 9 ) 

The difference in pressure between the two sides of the body is 
from equation ( 8 ) 


Ap = pVQ^Vg - (-PV1V2) = 2pv^Vg (10) 

and the total side force per unit length dF^^ds is equal to the average 
pressure difference ^ times the wetted hei^t h or 

dF — — 

= Aph = 2phvj^V2 (11) 

In order to evaluate the quantities v^ and V2^ consider first the 

limiting case of a body of infinitesimal transverse slope (e-> 0 ). For 
this condition the body is substantially reduced to a flat plate (see 
i*ig» 3) so that the flow in the horizontal direction (fig. 3 ("b)) is 
substantially undisturbed or 


V2 = (12) 

For a flat plate which is completely submerged the velocity distribution 
on the surface of the plate due solely to the vertical component of flow 
is (see, for example, ref. l) 
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VI 



T]/c 



(13) 


According to Wagner (refs. 2 or 3) the flow on a non- chine- immersed 
body of infinitesimal slope is substantially the same as that for the 
completely submerged flat plate given by equation (l3)* (The velocity 
variation for the case where the chines are immersed and where the upper 
surface is not wetted is considered in a later section of this paper.) 
Substituting equations (l2) and (13) into equation (8) gives the pres- 
sure due to yaw as 


P - Pi = IpV^Vtj 



(14) 


and the yaw or side force per unit length is then obtained from equa- 
tions (11), (12), and (13) as 


ds 


= 2pV^VT)h 






Evaluation of the integral in equation (l5) gives 


(15) 


ds 


2pV^V,^h 


(16) 


The relation between the wetted hei^t h and the penetration of the 
body 5 (see fig. l) has been obtained by Wagner for the unyawed case 

by using the flat-plate approximation. The resulting equations are 
given in references 2 and 3- For the special case of a straight-sided 
wedge Wagner obtained the relation 


h _ jr 

I"2 


(IT) 
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These same relations* are assumed to hold for the yawed case. Equations 
were also derived by Wagner for the sjnnmetrical pressure distribution Pj_. 

(See ref. 1 , 2, or 3 -) For the special case of a stral^t-sided wedge, 
the chines of which are not covered by water, Wagner’s equations reduce 
to the relation 



First approximation for the straight- sided wedge .- A first approxi- 
mation for the pressure distribution on a yawed strai^t-sided wedge is 
obtained by combining equations (l 4 ) and (18) as 



(19) 

The wetted semiwidth c which is equal to h cot P for a straight- sided 
wedge (see fig. 4 ) is then given by equation (17) as 

c = h cot P = ^ ^ cot P (20) 

The side force per unit length is given by combining equations (16) and 
(17) as 


dFri 

ds 

The center of pressure of the side force 


•cp 


/ 

Jo 

f 


Ap T) dT| 


Ap dT| 


( 21 ) 

is given by the relation 

( 22 ) 
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Substituting equations (lO), (l2), and (13) into equation (22) gives 


or 


'Hep “ 


T]/c 

I 2pV,V^ - pi' ■ , n dTi 

° _jnw 




jt/c 




dn 


(23) 


icp 

c 


10 




=rdff 

2 A'=' 


f- 


JlZs 




(24) 


and integrating equation (24) gives 

l££ = £ 

c 4 

The rolling moment per unit length about an axis at distance 


(25) 


a I above 


the keel (see fig. 4) is obtained as follows: The rolling moment due 

to the force on the right side of the wedge is equal to 

Fright sin S - 

where is the total force per unit length on the right side (which 

must be normal to the surface of ,the wedge) and where a^ sin P - 

is the center-of-pressure distance (see fig. 4). Similarly, on the left 
side of the wedge the rolling moment is 
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and the total rolling moment per unit length is 


right 


- Fieft) (a I sin p 


_ncp_'\ 
cos P/ 


or 

(Wt » - ^i*ft “1” 3 ) (n - 

■^1 

The side force per unit length is, by definition, equal to 

ds 

Fri^t sin P - Fpeft sin P; therefore the ro2_ling moment per unit length 
becomes 


dM| 

ds 



^cp \ 
sin P cos p/ 


(26) 


and substitution of equations (20), (2l), and (25) into equation (26) 
gives 


dM| 

ds 





A 

8 sin^p j 


(27) 


Second approximation for the strai^t- sided vedge .- This section 
is concerned with the development of a better approximation for the 
press\ire distribution and yawing forces and moments on a straight- sided 
wedge than is given by equations (l9)^ (2l) , (25), and (27) which were 
obtained for the limiting condition of a wedge of infinitely small slope 
(p->0). For this approximation, the unyawed pressure p^^ is obtained 

from equation (18) in the same manner as for the first approximation so 
that the pressure distribution may be expressed by combining equations (8) 
and (18) as 



■St cot p 1 




±2 


Za li 


(28) 
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The iDeBt method of estimating the quantity ’ which is the 

velocity distribution on the unyawed wedge, appeared to be the use of 
the detailed graphical solutions of Pierson for the constant-velocity 
impact of a strai^t-sided wedge on a smooth water surface. (See ref. 4 .) 
The resulting curves for which can be obtained from figure I6 

of reference 4 , are shown in figure 5 of this paper. It is noted that 
these theoretical curves are easily interpolated. Also shown in this 
figure is the theoretical variation of v^^V^ predicted by equation (l3) 

for the first approximation. 

In order to determine the quantity 'the solution is used 

for the symmetrical flow about a completely submerged straight- sided 
wedge. This problem has been solved by R. J. Monaghan (not generally 
available) and the method of solution is indicated in reference 5 * The 
resulting equation for is 



1 - 


(29) 


where \ is defined by the relation 


1 P 



and where F designates the gamma function. _The variation of Vg^'V-^ 
with r\/c as computed from this equation is shown in figure 6 for 
various angles of dead rise. Also shown is the variation predicted by 
equation (l2) for the first approximation. 


Curves of the product were obtained by multiplying 

together the corresponding values of i'rom figures 5 

and 6 and are shown in figure 7 together with. the prediction of the first 
approximation. These curves, which represent the pressure due to yaw 
(see eq. (8)), can be used with equation (28) to predict the pressiire 
distribution on the yawed wedge. 

The curves of figure 7 were graphically integrated to obtain the 
average value vjv^ and the results were substituted into equations (ll) 
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and (17) to ottain the side force per unit length. The resulting equa- 
tion for the side force is 


J = BiorpV^ (30) 

where the quantity B]_, which represents the ratio between the side 
force for the second approximation and that for the first approximation, 
is shown in figure 8. 

The center of pressure of the yaw force is found by graphical • 
integration of the pressure according to equations (lO) and (22) with 
the aid of the cirrves of figure 7- The resulting center-of -pressure 
locations are given by the relation 

(31) 

where which represents the ratio between the center-of -pressure 

distance for the second approximation and that for the first approxi- 
mation, is shown in figure 8. 

The rolling moment per unit length according to equations (20) , 

(26), (30), and (31) is 


^ = Bi.pVjT.sL - 


8 sin^P / 


(32) 


Peak pressure on the straight-sided wedge .- According to Wagner 
(ref. 3) the peak pressure on a straight- sided unyawed wedge of small 
dead-rise angle is approximately equal to the dynamic pressure corre- 
sponding to the side velocity of the edge of the wetted semiwidth or 



^ = Yt- cot P 
dt ^2 


where 
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For the yawed case the side velocity of the edge of -the wetted semiwidth 
is equal to the rate of expansion of the wetted width dc/dt plus or 
minus the side velocity of the body Vjj so that the peak pressure for 

this case becomes approximately 

and when dc/dt is assimed to be the same for both yawed and unyawed 
cases 


Pp = |p(v5f p * 

' cot e i ^ (33) 

Effects of chine Immersion for bodies of arbitrary shape .- In this, 
section consideration is made of the effects of chine immersion on the 
side force and rolling moment during a yawed impact. According to. equa-. 
tions (11) and (12) the side force is given approximately by the relatioir 


dF 


or in a more convenient form 


dF. 

- 2B2PhV^V5 


(34) 


where 



For the non-chine-iramersed condition Eg can be obtained as a first 

approximation by taking the average value of equation ( I3) . For the 
chine-immersed case, however, eqiiatlon (l3) does not apply and some 
other means must be used to determine This determination may be 

made approximately as follows: For deep immersions the fluid flow about 
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aji impacting body is very similar to the steady, separated, free- streamline 
flow about the same body. (See fig. 9») Solutions are available in 
reference 6 for the velocity during such fhee-streamline flow which indi- 
cate that B 2 is equal to approximately 0.3^ 0.^, and O .5 for dead-rise • 

angles of 0°, 20°, and 40°, respectively. 

The center of pressure (see eqs. 10, 12, and 22) according to the 
velocity relations of reference 6 is given by the relation 

T]p = i Egc = ^ E2b where E2 as 0.8 for dead-rise angles up to 40°. The 

rolling moment is obtained from this relation and equations (26) and (34) 
and becomes 

dMt / jtbEg \ 

If = g ,,3 (35) 


ANALYSIS OF OBLIQUE YAWED SEAPLANE LANDINGS AND PLANING 
FOR RESTRICTED SIDE MOTION 
General Analysis 


> The problem to be treated in this section is the fixed-trim, oblique, 
smooth-water landing of a seaplane of arbitrary constant hull cross 
section having a fixed angle of yaw. The angle of yaw if is defined as 
the angle between the landing plane (XZ-plane of fig. lO.(a)) and the 
plane of symmetry of the seaplane (X'Z' -plane of fig. 10(b)). The motion 
of the seaplane is assumed to be restricted to the XZ-plane; that is, 
there is assumed to be no side motion in the Y-direction. In an actual 
seaplane landing the seaplane is, of course, not restricted to motion in 
the XZ-plane but is free to accelerate perpendicular to the plane of 
landing (in the Y-direction) as a result of the side force. However, 
the restricted case corresponds to the model testing conditions for the 
experimental data of this paper aud it happens that the equations for 
the unrestricted- side-motion condition can be conveniently expressed in 
terms of those for the restricted case; therefore, it is expedient to 
consider the restricted case first. The unrestricted or free case is 
considered separately in a subsequent section. 

Consider the motion of the seaplane at a given instant of time (see 
fig. 10 (a)). The yawed float moves forward (in the X-direction) at a 
velocity x, vertically down (in the Z-direction) at a velocity z, and 
to the side (in the Y-direction) with zero velocity. If this motion is 
viewed perpendicular to the plane of symmetry of the seaplane (see 
fig. 10 (b)), it appears that an unyawed seaplane is moving forward (in 
the X' -direction) with a velocity x cos i)r, down (in the Z' -direction) 
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at a velocity z, and sidewise (in the Y’ -direction) at a velocity 
where 

Vt) = X sin Tif (36) 

In the theoretical treatment of the landing of unyawed seaplanes 
(refs. 7 to 11), it^s generally assumed that the flow processes in a 
given transverse plane, stationary in space and oriented normal to the 
keel (see aa' and hh' in fig. 10(b)), can be considered to be similar 
to those in the corresponding two-dimensional .problem. In the present 
case the velocity of penetration into this plane (aa' or bb') or the 
velocity normal to the keel is seen from figure 10(b) to be eqiial 

to 


= z cos T + X sin t cos ilf (37) 

and can be satisfactorily approximated by the simpler relation 

» z cos T + X sin t = ^ ( 38) 

for most practical angles of yaw since the cosine of the yaw angle is 
usually close to unity. The side velocity (parallel to the Y'-axis) is 
given by equation (36). 

EqiJiations (36) and (38) used together with the equations of the 
previous section on the two-dimensional problem provide a first approxi- 
mation for the side forces in terms of the instantaneous velocities i 
and . z and the normal displacement ^ . In order to determine these 

quantities as a function of time, the motion of a yawed seaplane in the 
XZ-plane is assiaaed to be the same as that for the unyawed seaplane; 
that is, the time histories of the draft, vertical velocity, horizontal 
velocity, vertical acceleration, and vertical hydrodynamic force are 
assumed to be the same for the yawed and xmyawed cases. (A subsequent 
section of this paper shows that this assumption is substantiated by 
experimental data for a float having an angle of dead rise of 22.5° for 
yaw angles at least up to 12°. It is also noted that this assumption 
is consistent with eq. (9)0 Procedures for computing the vertical 
motions for the unyawed case are given in references 7 to 12. Refer- 
ences 8 and 12 consider the seaplane of arbitrary constant cross section 
references 11 and 12 consider the V-bottom seaplane with considerable 
chine immersion, and references 9, and 10 deal particularly with the 
V-bottom seaplane having little or no chine immersion. 
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The Wide V-Bottom Seaplane 

The special problem of the yawed impact of a wide, constant-dead- 
rise V-bottom surface, the chines of which do not penetrate the water 
surface, is now considered in detail. In addition to the assumptions 
of the preceding sections the weight of the seaplane is assumed to be 
counterbalanced by an equal constant wing-lift force throughout any 
landing. The motions of the vmyawed seaplane for this case are dis- 
cussed in detail in reference 9 and equations and dimensionless plots 
for these motions are presented therein. 

Side force .- When equations (30), ( 36 ), and ( 38 ) are combined, 
the side force per unit length of float is obtained as 

= Biirp^x^ sin Tjr (39) 

ds 

The relation between ^ and s (see fig. 10(c)) is given by the 
relation 

5 = s tan T (4o) 

Combining equations (39) and (4o) and integrating the result gives the 
total side force as 


F 





Bj^itp^xs sin ilf tan t ds 


• • _ 

Bi^rp^xz^sin i|f 


2 sin T cos T 


(41) 


Since equation (4l) is based upon many inexact assumptions, it is unlikely 
to be always in good agreement with the experimental data. Consequently, 
equation (4l) is modified empirically here by the substitution of some 
empirical quantity B for Bj^ into equation (4l) so that 


Birp^xz^sin i|f 
2 sin T cos T 


(42) 


The empirical quantity B is probably substantially a function of dead 
rise alone and of the same order of magnitude as B]_. 
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In order to convert‘~the side force Into a dimensionless variable 
similar to those used in reference 9 > equation {k2) may first be 
rewritten in the form 


X 


BrtpXg^Q sin T)r 


6^ 


2/3 



[f(P)l^qp(A)itpg i 
6w sin T cos^T 


2/3 


2^tan T 


^9(a) xpgj 


(i^3) 


In the terminology of reference 9 the term 


[f(p )3 ^cp(A)rt'p^ 
6w sin T cos^ 


is the square of the dimensionless draft coefficient Cd.. The term on 

the left of the equation is defined here as the dimensionless side-force 
coefficient 


Cn - 


BitpXo^o sin t 




^ 2/: 
■> 


2^tan T |jf(3)]2(p(A)_pJTgJ 


. m 


so that equation (43) can be rewritten as ‘ 

c (45) 

' ^o Co - 

• • 

According to the analysis of reference 9 the time histories of C/Cq 
and C( 3 _ are functions only of the approach parameter k where 

ic = cos( 7 q + t) . The time history of xAo^ although easily 

sin 7o 

computed, is not a fimctlon of k alone. However, during many impacts 
the change in horizontal velocity is small; therefore, it appears 
desirable, in the interest of simplicity, to ignore this variation and 
approximate equation (45) by the simpler relation 
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so that the time history of the side-force coefficient becomes a func- 
tion of K alone. (This approximation leads to a conservative estimate 
of the side force.) Theoretical time histories of Ct| (according to 

eq. (^)) for various values of k have been computed from the solutions 
of reference 9 snd are shown in figure 11. The maximum values of 

are plotted against k in figure 12. This latter theoretical curve 
can be approximated within about 3 percent by the following empirical 
equation: 


C 


■Hmax 


10 

19 + 29 k + 


(0 < K < 10) 


The time to reach maximum side force is shown in figure 13 together 
with the times to reach TnAviTmim vertical load and maximum draft. 

Rolllpg moment .- The rolling moment can be obtained by substituting 

equations (38), (38), and (i«D) into equation (32) and integrating the 
result over the float bottom. Also B is .substituted for The 

res’olting equation for the rolling moment is 


= Bjtp^x sin if / 

Jo 


:/ sin T / tan t 

s tan T (at - s| ds 

8 sin2p 


Bnp^x sin if f&^z^ 


3in T cos T\ 2 24 sin^3 cos tJ 


Brtp^xz'^ sin if 


jt?Exz 


f (^ “ Q I 

2 sin T cos T \ 12a ^ sin^P cos tJ 


(47) 


and. substituting equations (42) and (44) into equation (47) gives 


M ^ — F^a. ^ (1 - 


rt^iz 

12a I sin^3 cos r) 


= C. 


a|Birp^QXo sin i 

2>^tan T 


Sii 


2/3 


jf(P)] 2(p(A)3tpg 


1 - 


rt^iz 

12a| sin^P cos tJ 
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The quantity Ei which represents the ratio between the center-of- 
pressure distance, for the second approximation ‘and that for ;bhe first 
approximation for the two-dimensional case is now replaced by the quan- 
tity E which represents the act-ual center of pressure. This quan- 
tity E is probably substantially a function of dead rise alone and 
is probably somewhat similar to Ej,. With this substitution and also 


by replacing z by f according to the defini- 

tion of C(i in reference 9 the preceding equation becomes 


where 


M| = Cti 


• • 

a^Bitp^^Xo sin t 




2/3, 




2 cftan T 


|f(3)] ^(p(A)rtpg 


1 - 


I^ECr 


sin T cos^T 


12a» sin^3 cos t l[fO)Pq>(A) «PgJ 


1/3 ' 


or 


Mp = 


jr%EW^Xo sin t 
4[f(3)]^9(A)g B±nh 


- Cd) 


m 


r. 


5g = a| 


1 2 sin^p 


[f(P)]VA)xpg 


6w tan T 


1/3 


(49) 


The rolling-moment dimensionless coefficient is now defined as 


'm 


^ sin t 

,4|f(P^^cp(A)g sin^P 


(50) 
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which from equation (48) is equal to 

Cux^ = - C(^ = - Ct|C(3_ (51) 

Since the time histories of both and are fimctions of k only 

but 5 1 is not a function of k, equation (51) shows that the time 

history of the rolling moment is different for each combination of k 
and 5^ so that a single family of generalized curves like figure 11 

cannot be given for the rolling-moment time history. For any specific 
case, however, the cii2rve of Cm| can be obtained by combining the 

pertinent Ct] curve of figure 11 (multiplied by 5|) with the corre- 
sponding Ct^C^ curve of figure l4- (which figure was obtained by com- 
bining the Ctj cvLTves of fig. 11 with the Cd. curves of fig. 6 of 

ref. 9) ■ The variation of the maximum rolling-moment coefficient with 
K is shown in figure 15 for several values of S^. 


Yawing moment .- The hydrodynamic yawing moment about ah axis per- 
pendicular to the keel at a distance a^ forward of the step ( see 

fig. 10(c)) is given by the relation 



ds 


a* - 


z 

sin T 


+ 


si ds 


(52) 


Substitution of equations {39) > (40), and (4-2) into equation (52), 
substitution of B for B]_, and integration of the result yields the 

following equation for the yawing moment: 


pz/sin T 

= Birpx^ sin tan t J s ^a^ - + s)ds 


Bnpx^z^ sin y / z 2z 

|a„ + 

2 sin T cos 


T V ^ sin T 3 sin t , 


= Fr^a Jl - 


3a^ sin t , 


(53) 
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In order to allow for the inaccuracies in the assumptions leading to 
equation (53)# this equation is modified here hy the insertion of an 
empirical constant G so that eqttation (53) "becomes 


Mf = FTia. (1 , 

^ ^ M 3a«. sin t 


(54) 


Substitution of equation (44) and the relation between 2 and 
into equation (54) gives the relation 


= 


BGWjcq^o ^ 

[f(3jJ^(p(A)g tan 


r [^(P)J ^'P(A)jtpg tan^T^ 


G 




Cd 


(55) 


When the dimensionless yawing- moment coefficient Cm^ is defined as 


and 5^ as 


ill 


BGWXq^q ,sin if 


5^ = 


[f(3)]2(p(A)g tan 


[f(P)]^<P(A)«Pg tan^rl 


V. 


6w 


equation (55) can be written in the form 


(56) 


(57) 


^m^ ~ ” C(i) — - ^T]^d 


(58) 


Since equations (51) and (58) for the rolling and yawing moments have 
been expressed in identical form, the time histories of the yawing mom- 
ents can be computed from the generalized cvirves of figures 11 and l4 
in the same manner as for the rolling moments. The maximum yawing mom- 
ents can be obtained from figure 15 in a like manner. 

Pressure distribution .- In order to make equations (l9) and (28) 
applicable to the three-dimensional problem, "V^ and Vtj in these 
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equations are replaced ty the values given by equations ( 36 ) and ( 38 ) > 
the term containing the acceleration is corrected by an aspect-ratio 
correction (see ref. 13) , and the angle of dead rise P is replaced 
by a more general quantity 6, called the effective angle of dead rise 
(see ref. l4) . Equations ( 19 ) and (28), respectively, become 



It cot 6 




2 £ctp(A) /7 

—r 



± 


X sin f c 



(59) 


for the first approximation and 



for the second approximation where 



(0 < A < “) (61) 


9 (A) 



(A >1.5) (62) 


and 


(Wetted length at keel)^ 

Wetted area projected normal to keel 


( 63 ) 


Reasons for the substitution of 6 for P, which substitution was 
Introduced by Pierson and Leshnover for the unyawed case in reference l4, 
are discussed in references l4 and 15» Several rather similar equations 
have been proposed for 6 in these papers and are as follows; 
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It cot &2_ 


It cot 01 


n cot ©2 


It cot 62 


K cot 03 


rt cot 03 


where 


= 2i 


/k ^ - 2 K sln^P - K^sln^ tan^r 
sin^3 + K^tan^T 




^3t cot 01 ^ 2^ 


iSh) 


..(^ zJ^p lL + cos2e 

sin23 + K^tan^T 


cot ©1 ^ 2^ 


= 2 . 


1 — ^ tan23 



^tan^T + ~ tan ^3 
^2 




(n cot ©2 ^ 2^ 


> ( 65 ) 


sin^t + — tan 23 cos^t 


cot ©2 < 2^ 


= 2 . 


1 - 


tan^T + 


sin^T + J^cos^T 


^ cot ©2 ^ 2^ 


> (66) 


(jt cot ©3 ^ 2^ 


J = - tan 3 

rt 


J = 0.293 


J = ^ tan 3 


(P->o) 




(3 = 22.50) ^ (67) 

(3 = 30 °) 


K R» 



3 tan ^3 cos 3 tan 3 sin^ 3 i 


1.7rt'= 


3.3Jt 


and 
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The relative merits of these equations for 6 are discussed in refer- 
ence 15. Equation ( 66 ) appeared to give the "best agreement with experi- 
mental data for the unyawed case for an angle of dead rise of 22 . 5 ° 

(ref. 15) and is therefore used subsequently for comparisons of theo- 
retical and e2q>erimental pressure distrihutions . 

Peak .pressTxre . - In oi^der to convert equation ( 33 ) for the peak 
pressure to the three-dimensional case, the terms V'q and are 

replaced hy x sin ilr and t,, respectively (see eqs. (36) and (38))j 
the following equation results: 


= /j cot 3 1 

V i j 


(68) 


Reference I6 shows that, for the \myawed case, equation ( 68 ) does not 
adequately take into accoimt the effect of trim on the peak pressure. 

In order to take this effect into account, equation ( 68 ) which, for the 
unyawed case, reduces to the relation 



(t = 0) (69) 


was empirically modified in reference 16 to obtain the following equa- 
tion for the peak~pressure at any finite angle of trim between 0° and 
90°: 



(t=o) (70) 


(See eq. ( 3 ) of ref. 16.) The corresponding modification of equa- 
tion {6Q) is 




(71) 


sin*=^T + 


COS^T 


'f cot p i 
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(Compare with eqs. (69) and (70)*) Reference 15 shows that equa- 
tion (70) was overconservative for the unyawed case for an angle of 
dead rise of 22.5° • This difficulty was remedied hy replacing the 

ir 1 

quantity — cot 3 "by an. empirical function where J is given by 

2 tj 

equations (67). Application of this same substitution to equation (7l) 
gives the following final equation for the peak pressure on a yawed 
wedge: 


Pp 

ipt 


sin^T + 


COS^T 



(72) ' 


Steady planing.- For the special case of steady planing (z = 0, 

^ = V sin r, x = V) the equations of the preceding sections for the 
side force and the rolling and yawing moments reduce to the following 
relations : 




BnpV^z^sin t 
2 cos T 

BjrpVSz^sin t 


a» - 


2 cos T 


Eit^z 


12 sin^p cos Ti 


(T3) 

(74) 




BrtpV^z^sln 
2 cos T 



Gz \ 
3 sin T/ 


(75) 


(See eqs. (4-2)^ (4-7), and (54), with E substituted for E;^ 

eq* (47)-) The ratio between the side force and the normal force is 
also of some Interest. According to reference 9 the normal force in 
planing is 


^ jt [f ( P)3 ^T>( A) pV^z^tan t 

h = 


2 


( 76 ) 
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so that when equation (73) is divided by equation (76) this ratio becomes 


Ft] _ B sin t 

[f'( 3 )l^q>(A) sin T 


( 77 ) 


The Yawed Seaplane After Chine Immersion 

Occurrence of chine immersion .- Chine immersion is considered to 
occur when the wave generated by the motion of the seaplane meets the 
chine-step intersection. From an examination of figure 10 (d) it is 
seen that this Immersion occurs when 


^ “ ^ch “ 


H cos T 
EL 


( 78 ) 


where R and L are wave-rise ratios illustrated in figure 10 (d) . 
For the special case of the strai^t- sided wedge ^ - \ 3 9^3. 
equation (78) becomes 


2ch 


b tan P cos t 
2 EL 


( 79 ) 


In the previous sections of this paper, the ratio R for the straight- 
sided wedge has been considered to be equal to jr/2 .(see eq. (17)) and 
the ratio L has been considered equal to unity so that equation ( 79 ) 
becomes 


b tan P COB t . 

zch = (80) 

jt 

Actually, according to experimental and theoretical information given 
in references 4 , 15 j 17 ^ and 18, the ratio R is either equal to 

or somewhat less than Jt/ 2 . The ratio L is either equal to or some- 
what greater than unity. Thus the two effects tend to compensate each 
other so that equation (80) may give a reasonable prediction of the 
draft at chine immersion. 

For cases of seaplanes which differ slightly from straight- sided 
V-shapes it is suggested that equation (80) be modified as follows: 
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^ch “ 


b tan Pi cos t 
n 


(81) 


where p]_ is the average angle of dead rise as indicated in figure 10(b) 
and is defined by the relation 


4. o 2H 

tan h = T- 


(82) 


Side force .- The side force after chine immersion on a seaplane of 
arbitrary consteint cross section is obtained by integrating equation (30) 
over the non-chine-immersed length of the seaplane (see eq. (1?))^ l>y 
integrating equation (3^) over the chine-immersed length of the sea- 
plane, and by substituting equations (36) and (38) into the resulting 
equation to obtain 


= SBpV^Vt, 


■'Zch/sin T 


h ds + 2B2 PV^Vt^ /• 


i/sln 1 

;h/sin 


h ds 


2pV^Vn(BSncl, + B2^nch) 

2px^ sin t + Ba^Sn^h) 


( 83 ) 


(The constant has been replaced by B in the first term of this 

equation to make tnls equation consistent with the analysis of the 
preceding sections for the non-chlne-immersed case at the instant of 
chine immersion.) The term is the wetted area of one side of 

the seaplane projected normal ,to the plane of symmetry of the seaplane 


at the time of chine immersion and AS- 


nch 


is the Increase in wetted 


area subsequent to chine immersion. These wetted areas include only 
the areas of the seaplane which lie below the plane of the chines. 
Usually the flow separates from the hull at the chines so that there 
is no wetted area above the chines. Cases where there may be such 
wetted areas, however, are considered in a subsequent section of this 
paper . 

The wetted areas needed for equation (83) can be obtained from an 
observation of figure 10(b) as 
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2H 

Tt tan T 


4it tan T 


and 


ASn^h =2^; 


2H 


sin T 1 C tan T 
b tan Pj./ z 


b tan 3i\ 


(tSin T TC tan t 


so that equation ( 83 ) becomes 


pb^x^ tan Pi sin tlr 

B tan Pi cos t 

sin T 

2ic 



tan 3i cos T 


1C 


)1 


(84) 


Rolling moment .- The rolling moment after chine immersion of a 

seaplane of arbitrary constant cross section, which depends on the 
characteristics of the individual hizll shape concerned, is not consid- 
ered in this paper. The rolling moment of a straight- sided wedge, 
however, can be easily obtained by integrating the rolling moment per 
unit length along the hull. For non-chine- immersed sections this quan- 
tity is given by equation ( 32 )j for chine-immersed sections, by equa- 
tion ( 35 ) • The resulting equation for the rolling moment is 


p^xb^tan P 

sin t 

B tan P cos 


icE 

2 sin T j 

It 

Vb 12 sin P cos 

2B h 


tan P cos T \ 

■ 

8 

sin P 

cos py\b 

” 1 





(85) 


(The constants and h^ve been replaced by B and E in order 

to make this equation consistent with the analysis for the non-chine- 
Immersed case at the Instant of chine immersion.) 
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Yawing moment .- The yawing moment after chine Immersion for a sea- 
plane of arbitrary constant cross section can he obtained similarly 
from integration of equation (52) as 




• ■ 

px^b^tan P]_ sin ilr 

B tan cos t 

!i . ° 1 

2 tan P2. 

sin T 

L 

b sin T 

^b 3 « 1 


L tan Pi cos t\ 

tk. 

G (z 

tan Pi cos 


K ) 

b 

2 sin 




( 86 ) 


Fiirther Refinements of the Chine-Immersion Theory 

The preceding development for the forces and moments on chine- 
immersed yawed bodies took into account only those forces acting on the 
parts of the bodies below the chines. On actual seaplanes, however, 
there are often substantially vertical walls extending above the chines 
which may also experience yawing forces and moments. In order to obtain 
a rough estimate of the order of magnitude of these forces, consider 
the problem of the steady planing of a rectangular flat plate of finite 
thickness w (fig. 16) . For this case the theory of the preceding 
parts of this paper gives no information regarding the imsymmetrical 
forces and moments. In figure l6(b) the motion of the plate is broken 
up into two components, one parallel to the plane of symmetry of the 
plate (i cos t) and one perpendicular to this plane (x sin ■>!<■) • Thus 
the side n]_ of the plate recedes from the water (is not wetted) and 

therefore has no water pressure on it. The side n2, on the other hand, 
may or may not be wetted. If it is wetted, it penetrates the water at 
the side velocity Vt) = x sin t and therefore has a water pressure on 

it. The side motion of this side n2 somewhat resembles the motion of 
a two-dimensional flat plate of beam w with separated flow behind the 
plate, for which case (see ref. 6 or 19) the average pressure per unit 
area is 


p = o.UpVn^ (87) 

strictly speaking, equation (87) should give correct results only for 
an infinitely large ratio of wetted length to plate thickness (two- 
dimensional flow) and for zero beam (such that there is no interference 
between the flow below the hull and on the sides) . If the ratio of 
wetted length to plate thickness is small, equation (87) for zero beam 
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may considerably underestimate the actual pressure and, if the beam is 
increased, the flow under the float will tend to separate at the sides 
and thus reduce the side pressure. Thus the two effects tend partly to 
compensate each other. With these reservations in mind, it is suggested 
that, in lieu of a more satisfactory procediire, a pressure equal to the 
value given by equation ( 87 ) be considered to act on one side of any 
parts of a seaplane above the chines which are submerged below the water 
surface and which are believed to be wetted by the water and that the 
yawing forces and moments on these submerged parts be obtained by inte- 
grating this pressure over the wetted areas concerned. Althou^ the 
resulting yawing forces and moments are hardly precisely correct, even 
for the case of the long rectangular flat plate, they are probably rea- 
sonable as a first approximation. 


Another Theory for the Side Force on Chine -Immersed or Non-Chine- 


Inmersed Straight-Sided Wedges 


The preceding sections of this paper have presented an analysis of 
the yawing forces based primarily upon the concept of two-dimensional 
flow in transverse planes. The applicability of the equations resulting 
from that analysis, however, has not yet been established for all prac- 
tical conditionsj therefore, it is desirable to consider also other 
possible methods for computing the yawing forces. For the special case 
of the side force only, such a method, which happens to be simpler in 
form than that of the preceding sections, is as follows. Consider the 
case of a planing wedge (fig. 17) • When this wedge is viewed in a plane 
perpendicular to the direction of motion, it is seen that the projected 
area of one side of the wetted surface (Sp in fig. 17) is larger than 


the projected area of the other side S 2 , the corresponding total wetted 
areas corresponding to Sp and Sp being assumed to be equal. If the 
total force (normal to the plating) on each side of the float (Fp and 
F 2 ) is considered to be divided approximately proportionally to this 


projected area difference 


r Fp 

that is ^ = 


the side force is equal 


F2 S2y 

to the difference between the horizontal components of these forces or 


sin p = |l(Sl - 

and the force normal to the keel is equal to the sum of the normal (to 
the keel) components of these forces or 


S2^ sin P (88) 


Fti = (Fp - F2) sin P = (Fp - 




F1S2) 

Si/ 
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= ^ (Si + S2) cos 3 (89) 

When equations (88) and (89) are combined, the side force can be 
rewritten as 


I 3 . = tan 3 (90) ' 

Si + S2 

By a consideration of the geometry of the float it can be shown that the 
Si - S2 

ratio is given by the relation 

Si + S2 


= (Fi + P 2) cos 3 = |Fi + 



cos 3 


^1 ~ ^2 ^ tan t tan 3 
Si + S2 sin T 


(91) 


which is valid for both chine- immersed and non- chine- im mersed wedges. 
Equation (90) therefore can be written as 


= tan t tan^3 
F^ sin T 


(92) 


Substituting the relations = x sin (see eq. (36)) and 

= X sin T cos Sr (see eq. (37)) and rearranging equation (92) yields 
the following relation: 


^ cot 3 = — tan 3 (93) 

P5 Vj 


The significance of -the right-hand side of equation (93) is indicated 
in figure I8 which shows transverse cross sections of a non-chine- 


Vr 


immersed wedge for values of ^ tan 3 less than and greater than unity. 
V ^ 

For -3. tan 3 larger than unity, only one side of the float is wetted 


by water; therefore, for these cases the ratio between the side and 
normal forces is equal to the tangent of the angle of dead rise (since 
the total force is substantially noimal to the plating) or. 
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cot P = 1 


(Yil tan 3 > 1 


( 94 ) 


Vti 

For values of — *• tan P less than unity, the data of figure 19 show 

that experimental planing data for a wedge having an angle of dead rise 

F 

of 22 . 5 ° give higher values of -2. cot 3 than are predicted by equa- 

F^ 

tlon (93)* (See the appendix for a description of these tests.) This 
result suggests that a more adequate equation for the side force could 
be written in the form 


^ cot 3 = f/^ tan 31 (95) 

VC / 

where the fiinctional relation is to be determined empirically. Since 
this relation must satisfy equation (94), it is probable that the rela- 
tion will not vary much with the angle of dead rise. The following 
simple equation is proposed for this relation (see fig. 19) 


— - cot 3 = sinj !- tan 3} 


F 







3^ 


( 96 ) 


Equations (94) and ( 96 ) , which were obtained for the planing case, 
may be considered to apply also approximately for the impact condition. 
This application is shown in figure 19 where experimental impact dat^ 
for a wedge having an angle of dead rise of 22.5° (see appendix) for 
various trims are shown to be approximately in agreement with the vari- 


ation of 


^ cot 3 


with — !• tan 3 given by the planing data and equa- 

tion ( 96 ) . Although it is realized that the actual variation probably 
changes with dead rise, trim, flight path, and other variables, it is 
believed that equation ( 96 ) will give reasonable results, at least for 
rough calculations, for arbitraiy yawed landing conditions. 


It is noted that for — !■ tan 3 near to and greater than unity the 

T?r» 

experimental impact data give values of 




cot 3 greater than unity. 
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This result indicates that there is a negative pressure on the low- 
pressure side of the float. Althou^ such an effect would not exist 
for the planing case ^the water would not wet one side of the float for 

during these particular yawed impacts, during earlier 

stages of the impacts (before the time represented hy the test points 

shown) the quantity ^ tan |3 was small enough to permit wetting of 

both sides of the float. The suction results from the tendency of the 



water to separate from the low-pressure side as the quantity 
Increases above unity. 


Vn 

^ tan P 


MALYSIS OF 0BLI5PE YAWED SEAPLABE LABDIBGS FOR FREE SIDE MOTION 


Basic considerations .- In the preceding sections of this paper it 

has been assumed that the motion of the seaplane during landing is 
restricted to the plane of its landing (XZ-plane of fig. 10(a)) or that 
there is no side motion (in the Y-direction% Since in an actual landing 
a seaplane is free to accelerate to the side as a result of the side 
force, it is desirable to determine the relations between the yawing 
forces in such a free yawed landing and those in the previously discussed 
restricted yawed landing. As for the restricted- landing analysis, it 
is assiimed that the vertical motions and forces are unaffected by the 
angle of yaw and that the change in horizontal velocity during an Impact 
is small. 


> Consider the motion of the seaplane in a horizontal direction 
normal to the keel. The force in this direction, which is the side 
force must (according to Newton’s second law) equal the mass of 


the seaplane 


H. 

g 


times the rate of change of the side velocity 


dVTlf 

dt 


or 


F 


Tlf - - 


W 

g dt 


(97) 


^the negative sign Indicates that Fj^^ and have been defined as 

positive in opposite directions^. The side force for the general case 
of a body of constant arbitrary hull cross section, whether restricted 
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in side motion or not, is given by equation (83) as 

r, = 2 pV5V^(bs„^^ + ®2^n^h) 

For the restricted case of the previous sections Vt| was considered 

to be substantially constant as a consequence of the assumption of no 
side motion and of the assumption i = thus was equal to its 

initial value or 

^Tla = ^Tlo ^98) 

and equation (83) for the restricted case becomes 

''Tla = + Bs&Snah) (99) 

For the free case V- is a variahle (having, of course, the same Initial 
value as for the restricted case^ and Ftj is given by equation (83) 

as 

Ft^^ = 2pV^VT,^(BSnch + Bg^ncn) (10°) 


When eqtiations ( 99 ) and (lOO) are combined, the relation between the 
side forces for the free and restricted cases becomes 


^^f ^Tlf 
^Tla ^^o 


( 101 ) 


^It should be noted that the variables and in 

eqs. (99) and (lOO), which depend only on the vertical and horizontal 
motions, have been considered identical for the two cases of free and 

restricted landings.^ A combination of equations ( 97 ) and (lOl) and 
a rearrangement gives the following relation (see eq. (38)): 
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VTlf ■" " 


^T]a 

Wxq sin t 


dt 


whicla, when integrated, becomes 



■/ 




Wxq sin t 


dt 


( 102 ) 


(103) 


When equations ( 103 ) and (lOl) are combined, the ratio between the 
side force for the free and restricted cases is 




'T)f 


= e 


'0 Wxq sin t 


dt 


(104) 


The same ratio is similarly applicable to the rolling and yawing moments, 
that ..is. 


\ ^ % 



_ ^1a 
Wxq sin T|f 


dt 


(105) 


Equation (l05) gives all the relations necessary for converting the 
equations of the previous sections for the restricted landing to those 
of the more realistic free landing. 

The wide non- chine -immersed wedge .- For the special case of the 
wide non- chine- immersed straight- sided wedge, equation (105) can (after 
much trigonometrical manipulation) be expressed in the following form: 







-a(l + k) 
e 



dCt 


(106) 
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where 


a 




[f(p)] 2 (p(A) 


(107) 


(see definitions of Ct^ and k) and where C^j is given hy fig- 

tire 11. (it should be noted that the quantity a is primarily a func- 
tion of the seaplane dead rise alone except for the quantity cp(A) 

which is usually near unity.) The integral (l + k) / dCt 

'^0 

(obtained from graphical integration of the curves in fig. ll) is shown 
in figure 20 , for various values of k, as a function of C-^. The 

reduction of maximum side force due to the change from restricted to 
free motion (obtained from fig. 20 and eq. (106)) is shown in figure 21 
as a function of k for various values of a. The maximum side-force 
coefficient for the free case, as defined by the relation 






Bjtpxpgo sin t 


6w 


2/3 




2 X tan T 


[f(^)] VA)npg, 


(108) 


(same definition as for the restricted case), is given in figure 22 as 
a function of k for various values of a. 


CONSTAMTS FOR THEORETICAL SOLOTIONS 


Dead-rise— aspect-ratio function .- In order to make theoretical 
solutions for the unsymmetrical forces and moments for the wide straight- 
sided wedge it is first necessary to solve the symmetrical problem 
according to the methods of reference 9 * This solution Involves the 

use of the quantity >^|f(p)pcp(A) which is related to the maximum 

vertical force (during a landing with wing lift equal to seaplane weight) 
in reference 9 by the relation 
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ho 




^ ^ sin t cos^t 

CZmax^o^V 


(109) 


where OzmBx i® given "by the solid ciirve in figure 23- This quantity 


'^[fO)] 2<p(a) 


was determined for an angle of dead rise of 22.5° from 


equation (IO9) hy using experimental maximum -load data from tahld I of 
this paper and from table II of reference 15 . The resulting values 
for all these impacts for which the maximum load occurred prior to 
chine I mm ersion are shown in figure 24. The solid faired curve drawn 
throu^ these test points was used for all computations in this paper. 


For angles of dead rise other than 22.5® it is suggested that the 
dead- rise— aspect-ratio fxmction be computed from the relation 


[f(3)]2<p(A) = - l)^(A) 


where qp(A) 


is given by equations (61) and (62) and where 


(no) 


A.S_5!£_§ (111) 

It tan T 

according to equations (17) and (63). Equation (llO) gives the dead- 
rise function proposed by Wagner in reference 2 modified by an aspect- 
ratio correction. Frcan a consideration of the experimental and theo- 
retical information in references 4, 5, 7, 9, 10, 15, and 20 and in 
this paper it appears that equation (llO) is reasonably satisfactory 
for angles of dead rise from 30® to 50®. For angles of dead rise 
smaller than 30® the experimental evidence in references "J, 15, and 20 
and in this paper indicates that the actual dead-rise— aspect- ratio func- 
tion is somewhat smaller than the value given by equation (llO) (see, 
for example, the dashed line in fig. 24) . The experimental evidence 
in these various references, however, is not completely consistent or 
conclusive; therefore, no attempt is made in this paper to define more 
accurately the dead-rise function. 

Yawing coefficients .- The following values of the quantities B, E, 
and G were empirically chosen from the experimental data in this paper 
and were used for all computations in this paper for p = 22.5°: 


B = 1.2; 


E = 0.7; 


G = 1.2 


( 112 ) 
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When it Is necessary to detenaine the quantity B for other angles 
of dead rise, it is noted that the value of B given hy equations (ll 2 ) 
is substantially equal to the value given' by the second two-dimensional 
approximation Bi shown in figure 8 . Consequently, it is recommended 
that this curve for B^ be used for B for computations for other 
angles of dead rise. 

It is noted that the value of E given in equations (ll 2 ) is 
smaller than the corresponding value of E]^ (see fig. 8 ). Purely 

arbitrarily it is suggested that the faired curve drawn throu^ 

E = 0.7 at P = 22.5 be used to compute E for other angles of 
dead rise. 

It is suggested that the value of G given by equations ( 112 ) 
be used for all angles of dead rise. 


COMPAEISONS BETWEEN THEORY AND EXPERIMENT FOR THE 
NON-CHIHE-IMMERSED WEDGE 


Experimental motions and forces as obtained from a yawed landing 
and planing investigation at the Langley impact basin (see appendix) 
are compared with the theoretical predictions of this paper in fig- 
ures 23 and 25 to 37 * In all cases the theoretical relations used are 
those derived in the section of this paper entitled "The Wide 'Y-bottom 
Seaplane . " 


Landing Investigation 

Vertical loads and motions .- Experimental time histories and maxi- 
mum values of the draft, vertical velocity, and vertical acceleration 
(or vertical hydrodynamic load) for landings with yaw angles between 
0 ° and 12 °, shown in figures 23 and 25 to 28, appear to be substantially 
Independent of the yaw angle and are in reasonable agreement with the 
predictions of the theory of reference 9 - 

Side load .- Comparisons of theoretical and experimental side loads 
shown in figures 29 and 30 indicate that the experimental side loads 
are in reasonable agreement with the theoretical values for yaw angles 
at least up to 12 °. 

Rolling and yawing moments .- Comparisons of theoretical and experi- 
mental maximum rolling moments, rolling-moment time histories, and 
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yawing-moment time histories, shown in figures 31^ 32, and 33 ^ respec- 
tively, indicate that these quantities are in reasonable agreement. 

Times to maxlpium vertical load, maximum side load, and maximum 
draft .- In figure 3 ^ comparisons axe shown between the experimental - 
and theoretical times bo maximum vertical load, maximum side load, and 
maximum draft. The experimental and theoretical times are seen to be 
in reasonable agreement. 

Pressure distribution .- Theoretical pressure distributions, as 
computed from the two approximations of equations (59) and (60), are 
shown in figure 35 together with the corresponding experimental pres- 
sure distributions for several trims and yaw angles. Fair agreement 
is seen to exist, the second approximation theory (eq. (60)) giving 
somewhat better results than the first approximation theory (eq. (59))- 

Peak pressures .- Experimental and theoretical peak pressures are 
compared in figure 36 for several yawed impacts . Fair agreement is 
seen to exist. 


Planing Investigation 

Experimental yawing forces and moments during planing are compared 
with the corresponding theoretical forces and moments (computed from 
eqs. ( 73 ) to ( 75 ) and ( 77 )) in figure 37 - It is seen that the experi- 
mental side loads, rolling and yawing moments, and ratios of side to 
normal load are in reasonable agreement with the corresponding theoret- 
ical quantities. 


CONCLUDING EEMAEKS 


An approximate theory has been developed for predicting the unsym- 
metrical forces and moments during the yawed impact or planing of sea- 
planes. From comparisons of the results of theoretical calculations 
with experimental side force, rolling and yawing moment, and pressure 
data presented in this paper for a non-chine- immersed straight-sided 
wedge having an angle of dead rise of 22.5° and yaw angles up to 12°, 
the proposed theory appears to give reasonable results for the cases 
tested. Although it is probable that the proposed theory is also at 
least qualitatively applicable to other angles of dead rise and to the 
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chine -immersed case, the extent of quantitative agreement for such 
conditions remains to he determined by future experiments. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 21, 1952. 
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APPENDIX 

EXPERIMENTAL INVESTIGATION 
Description of Tests and Instruments 


A series of smooth-water yawed and unyawed landing and planing 
tests was made at the Langley impact basin (ref. 21) with a prismatic 
V-bottom float having an angle of dead rise of 22.5°. Details of this 
model are shown in figure 38 * 

Landing tests .- Thirty-six landings were made with the float loaded 
to a weight of 1177 pounds which corresponded to a beam- loading coeffi- 
cient of 0.^. These landings covered yaw angles between 0° and 12° 
and initial flight-path angles between approximately 3° and 9°* During 
each landing a compressed-air engine (described in ref. 21) exerted a 
vertical lift force on the model equal to its wei^t so that the model 
simulated a seaplane with wing lift equal to the weight of the seaplane. 
Otherwise the model was free to move in the vertical direction. The 
model was attached to a towing carriage weiring approximately 5^^ pounds. 
Because of this large additional carriage inertia, the model did not 
slow down very much (horizontally) during any landing. 

Time histories of the following quantities were measured during 
each landing: vertical draft, vertical velocity, vertical acceleration 

(vertical force), horizontal velocity, side force, rolling moment, yawing 
moment, and pressure distribution. The instruments used for measuring 
horizontal velocity and initial- vertical velocity are described in ref- 
erence 21. The instruments used for measuring time histories of the 
vertical motions and of the pressure distribution are described in ref- 
erence 15 . The locations of the pressure gages used in this investiga- 
tion are shown in table II and figure 39 * The hydrodynamic forces due 
to yaw were measured as follows: A six-component strain- gage -type dyna- 

mometer measured the side force, rolling moment, and yawing moment at a 
position above the float. These moments and forces differ from the 
corresponding hydrodynamic quantities because of inertia forces and 
moments resulting from structural oscillations of the float beneath the 
dynamometer. The inertia force due to side oscillations, which occurred 
at a frequency of approximately 3-5 cycles per second, was measured by 
a side accelerometer mounted on the float. The effects of another struc- 
tural oscillation, which occurred at a frequency of approximately 13 cycles 
per second and which can be seen in figmre 33 ^ were hot taken into account. 

Planing tests .- Four planing runs were made in the Langley impact 
basin according to the procedure outlined in reference 22. During these 
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runs, measurements of tlie pressure distribution, the side load, the 
rolling moment, the yawing moment, the total force normal to the keel, 
and the horizontal velocity were recorded. The instruments used were 
the same as for the landing tests except that the normal force was 
measured by the previously mentioned dynamometer for the planing runs. 
The draft corresponding to each planing run was obtained by observing 
the wetted length, as given by observation of the experimental pressiire 
distributions, and multiplying it by the sine of the trim. If the water 
rises up in front of the float, this quantity (wetted length times the 
sine of the trim) may be greater than the actual draft; however, this 
rise is probably small for the conditions tested. 


Precision of Measurements 

The instrumentation used in these tests gives measurements that 
are estimated to be accurate usually within the following limits: 


Horizontal velocity, ft/sec iO.5 

Vertical velocity at water contact, ft/ sec ±0.2 

Vertical velocity after water contact, ft/sec ±0.5 

Vertical acceleration, percent ±5 

Pressure, Ib/sq in ±2 ±0.1p 

Time, sec ±0.005 

Draft for landings, ft ±0.03 

Side force, lb ±50 

Rolling moment, ft-lb ±100 

Yawing moment, ft-lb ±200 


In order to evaluate properly the experimental data from most 
instruments, the dynamic response characteristics of the instruments 
and the .corresponding recording galvanometers had to be taken into 
account. Analysis of this dynamic response for the experimental condi- 
tions of these tests showed that the response characteristics of the 
instrument-galvanometer circuits were of such a nature that the magni- 
tudes of the recorded motions and forces should be reasonably accurate 
but that the circuits had small time lags. All the data presented in 
this paper have been corrected for these time lags. 


Experimental Results 

Landing tests .- The initial vertical, horizontal, and resultant 
velocities,- fli^t-path angles, trims, and yaw angles for all landings 
are presented in table I together with the measvired over -all loads 
and motions. Time histories of these quantities for several landings 
are shown in figures 26 to 28, 30, 32, and 33* The values of the maxi- 
mum pressures recorded on each pressure gage are presented for all runs 
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in table III and the corresponding times of maximum pressure are pre- 
sented in table IV. Time histories of the vertical velocity are pre- 
sented in table V. In table VI instantaneous pressure distributions 
from several landings together with the corresponding measurements of 
time, draft, vertical velocity, and vertical load factor are given. 
The corresponding horizontal velocities are substantially the same as 
the initial values given in table I since the change in horizontal 
velocity during any impact was small. 

Pinning tests .- The experimental over-all loads and pressure- 
distribution measurements obtained during the four planing runs are 
presented in table VII together with the planing conditions of trim, 
yaw, draft, and horizontal velocity. 
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TABLE VI - Concluded 


INSTAWTAKEOUS HffiSSURE DISTRLBUTIONS - Concluded 



Pressure (Ib/sq in.) at gage nuiBber - 
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Figure 1 .- Ttie oblique in5>act of a symmetrical body. Two-dimensional problem. 



(a) Vertical conqjonent, (b) Horizontal component. 

Figure 2 .- Resolution of the oblique impact into two component flows. 

Two-dimensional problem. 



MCA TN 2817 


61 





♦ 





(a) Vertical component, (b) Horizontal component. 

Figure 3,- Resolution of the oblique impact into two component flows at 
very ama.l 1 body slopes. Two-dimensional problem. 


c 



Figure i)-,- Geometric relations for a straigb't-sided wedge. 


p 

(deg) 



Figure 5.- Variation of the velocity ratio on the side of a 

straight- sided two-dimensional wedge. 
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(c) Yawed landing of a non-chine -immersed seaplane without longitudinal 
water pile-up, viewed perpendicular to seaplane plane of symmetry. 



(d) Yawed landing of a non-chlne-immersed seaplane with longitudinal 
water pile-up, viewed perpendicular to seaplane plane of syrmetry. 


Figure 10 .- Concluded. 




0 .5 1.0 1.5 2-0 2.5 3.0 3.5 U.0 U.S 5.0 


(1 +k) Ct 

Figure 11.- Theoretical time histories of the side-force coefficient 
for a non- chine -immersed wedge. (Ordinates and ahscissas have “been 
multiplied by the factor 1 + k to facilitate interpolation. ) 
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Figure 12 .- Theoretical variation of the maximum side-force coefficient 
with the approach parameter for a non-chine-immersed wedge . 
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0 1 2 3 U 5 6 7 8.9 10 

Approach paraaetar, k 

Figure 13 ,- Theoretical variation of the tijne coefficient with the approach 
parameter for a non-chlne-immersed wedge. 
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(a) Curves for small values of or 6^. 


Figure 15.- Theoretical 
coefficient with the 


variation of the maximum 
approach parameter for a 


rolling- or yawing-monent 
non-chine -immersed wedge. 
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Approach parameter, k 



(b) Cutves for large values of 5 ^ or 6 ^ 


Figure I5.- Concluded. 
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(a) Rectangular flat plate viewed perpendicular to the plane of motion. 



(b) Rectangular flat plate viewed perpendicular to the plane of symmet^ 

of the plate. 


Figure 16.- The yawed planing of a rectangular flat plate. 
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Side vlcTTr 

Flgiire 17 •- Yawed wedge during steady planing. 




Figure 18.- Transverse ' sections of a non-chlne-lnimersed wedge for values 

Vt, 



Figure 19.- Variation of ^ cot p with ~ tan p for a straight -sided wedge, 
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’sp \ >1 ♦ T) 



r>Ct 

Figure 20 .- Theoretical variation of the quantity (l + k) / C« djCt 

' ^0 

vitli the tiaie coefrident for variouu values of approach paraineter 
for a non-chine -Immersed vedge.' (The abscissa has been multiplied 
by 1 + K to permit use of this curve with figs. 11 and lit.) 
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IfaxliftiuB Tcrtical-load-f actor coefflolent 



Approach paraiiatert k 

Figure 23.- Comparison teiveen tlieoretlcal and experimental variations 
of tbe maximum vertical "load-factor coefficient wltli the approach 

^CU. CUUC l/CJ. • 
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0 .Oil .08 .12 ,16 .20 

Time, t, sec 


(a) T = 3 . 2 °j Zq = 4.2 feet per secondj Aq = 71.5 feet per second. 



0 .Oil .08 ,12 .16 .20 


Time, t, sec 

(t) T = 9*3°i Zo = 3-3 feet per secondj Xq = ^T.l feet per second. 

Figure 27.- Comparison between theoretical and experimental vertical- 
velocity time histories, p = 22 . 5 °j W = 1177 pounds. 
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Theoretical solutions 

(deg) 

O Experimental data 


va 


Time, t, sec 

k ,2 feet per second: xo = 71-5 feet per second 


Time, t, sec 

(b) T = 9 - 3 °J zq = 3.3 feet per second) Xq = 67. 1 feet per second. 

Figure 28 .- Comparison between theoretical and experimental vertical 

load time histories. P = 22. 5° J W = 1177 pounds. 





side -force coefficient 


% 
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Theoretical solution from figure 12 

(deg) 

O 3] 

□ 6 I Experimental data from table I 

2 P - 22.5°; T - 3.2°, 6.3°, 9.3° 

C5^ c5 After chine immersion 


Approach parameter, k 

Figure 29 .- Comparison of theoretical and experimental variations of 
maximum side-force coefficient with approach parameter. 





^n^^^nfy.Illolllent coefficient 

Rolling-fflosisnt parameter 



Flaure Comarlsott of theoretical and experiiKnt^ 
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2 .3 

Time, t, sec 


(b) T = 9.3°j Zq = 3-3 feet per secondj .1 feet per second. 


Figure 32.- Comparison between theoretical and experimental rolling -moment 
time histories. 3 = 22.5°j t = 12°} a| = 2.96 feetj ¥ = II77 pounds. 




Tina coefficient 
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2.1» r- 


2.0 L 


Theoretical solutions from figure 13 

O At max. vertical load] Experimental data from table I 


□ At max. side force r 

O At max. draft J 

o' After chine immersion 


p - 22.5°: T - 3.2°, 6.3°, 9.1 
y - 3°, 6^, 9 °, 12® 


1.2 U 


At max:, draft 
At max. side force 
- At max. vertical load 




Approach parametw, k 

Figure 34 .- Comparison of theoretical and experimental time coefficients 

for a non-chine -immersed wedge. 
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c 20 


— Equation ( 72 ) (Vertical velocity from table V) 
O High-pressure side! Experimental data from 
□ Loir-pressure side tables III and IV 


High-pressure side 


tables III and IV 


Loir-pressure side 


(a) Run 5 : T = 3 ‘. 2 °j f = 


.08 ,12 .16 

Time, t, sec 

12 ° j Xq = 71*5 feet per second j Zq 
per second. 


= 4.2 feet 


cr 20 


V 

§ 10 



High-pressure side 


Loir-pressure side 


Time, t, sec 

(b) Run 27: T = 9 - 3 °J = 9 °J Xq = 67.2 feet per second; 

Zq = 4.3 feet per second. 

Figure 36.- Comparisons of computed and experimental peak pressures. 





VO 

ON 




(a) Side force, Fvj, 


(b) Rolling monjeiit, M|, 


poimds. 


poimd-feeb. 




(c) Yavlng moment, M^, (d) Ratio of side to noniial force, 

pound-feet. 




Figure 37* - Comparison of experimentel ari<^ theoretical forces and moments 

during yaved planing, p = 22.5'^. 


HACA TN 2817 





HACA-LloslaT - U-U<Kt - ^ WQ 


VD 

00 


1 



t 



oU 

08 

0 

H 









02 U 

* 


03 

07 

ul2 

017 

023 






016 

022 



02 

06 

oil 

021 

026 

028 


90 0ii. 

®2u 

025 


ol 

05 

010 

015 

019 

027 




018 




033 035 
029 036 

nn 

j\y 

031 037 

032 031 ^ 


f 

> 


1 

O ^-in»-“dlameter gages 
* p-in,- diameter gages 



Chine 



/ 

^ Keel 



Chine 





: c . • 




^ I 


d A 


NACA TN 2817 



